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PREFACE 

Multispectral scanners onboard NASA unmanned Landsat satellites 
provide an ideal source o£ current data for earth resources 
applications. The Detection And Mapping (DAM) package was origi- 
nally developed at the Johnson Space Center for rapid conversion 
of the Landsat digital data into hydrographic maps matching 
standard topographic quadrangle series. Recent improvements in 
both the manual procedures and computer programs within the DAM 
package make it easier to use, faster, and more general purpose. 

Documentation and software for the DAM package are available to 
all public and private agencies, in accordance with the NASA 
policy of encouraging maximum use of remote sensing technology, 

Published documentation, of which this is volume 3, is comprised 
of the following volumes: 

Volume 1: General Procedure 

Volume 2: Software User Manual Cin two parts) 

Volume 3: Control Network Establishment 

These volumes supersede the previous documentation published in 
1973. Software releases prior to version 7602 cannot be used 
with the current documentation. 
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1. INTRODUCTION 

The Landsat satellites circle the earth every 103 minutes in 
near-polar, sun- synchronous orbits. During each southbound 
(daylight) pass the multispectral scanner (MSS) onboard can pro- 
duce a digital image of the earth's surface radiance in a con- 
tinuous swath 185 kilometers (100 nautical miles) wide. 

The MSS digital data are recorded in a dynamic, scanner-oriented, 
coordinate system. To allow mapping and direct correlation with 
information from other sources, MSS data must be transformed into 
a static, earth-oriented coordinate system. The approximate 
coefficients for such a transformation may be derived from system 
measurements of spacecraft orbit and attitude and of scanner geom- 
etry. Exact transformation coefficients must be computed from a 
network of control points for which both scanner-oriented and 
earth-oriented coordinates are known. 


The remainder of this section describes the MSS data characteris- 
tics and the criteria for establishing a control network. 


1.1 MSS CHARACTERISTICS 


The MSS records spectral information .n four channels covering 
the following wavelengths; 


Channel 

1 

2 

3 

4 


Spectral Band 

4 

5 

6 
7 


Wavelength (micrometers) 


0.5 - 0.6 j 
0.6 - 0.7 \ 
0.7 - 0.8 ) 
0.8 - 1.1 


visible 

reflective 

infrared 


The MSS uses an oscillating mirror to scan perpendicular to the 
satellite's orbital path. The instrument's instantaneous field 
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o£ view is sampled 3,240 times* during each active scan line and 
the analog detector signals for each of these picture elements 
(pixels) are digitized. Satellite motion provides along-track 
progression of the scan lines. Rotation of the earth under the 
satellite orbit results in progressive offset to the west of 
successive scan line coverage. 

1.2 MSS DATA PRODUCTS 

The continuous MSS coverage is t ansmitted back to the ground and 
segmented into scenes of 2,340 scan lines by 3,240 samples* 
(approximately 100 by 100 nautical miles on the ground) or 
7,581,600 pixels. Each scene is then further processed into a 
number of data products available to the public. The two products 
which concern us in this volume are; 

• System- Corrected Images . Each image covers an entire scene, 
except for 42 lines at the top and 42 lines at the bottom. 

The images are radiometrically calibrated and approximately 
corrected for scanner geometry, spacecraft altitude/attitude, 
and earth rotation based on satellite tracking and telemetry 
data. Figure 1 shows a typical system-corrected image of 
channel 4 (band 7). The ERTS-1 Data User's Handbook (ref- 
erence 1) describes the various scales and formats available. 

• "System-corrected" Computer-compatible Tape (CCT ) . Digital 
MSS data for each 100 by 100 nautical mile Landsat scene are 
divided into four 25-mile wide-strips (numbered 1 through 4 
from left to right) . The data are radiometrically calibrated, 
but not corrected for scanner geometry or orbital characteris- 
tics. The data for all four strips of the scene are then 
recorded on either 1, 2, or 4 reels of computer tape, depending 
on the file structure and recording density. 

*This is an approximate number. The number of samples per scan 
line is constant within any individual scene, but may vary from 
scene to scene. 
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1 . 3 SCANNER- ORIENTED COORDINATES 

Based on the conventions explained above and in reference 1, the 
following standardized notation will be used for identifying the 
scanner-oriented coordinates of a pixel: 

Cl) Scan lines are numbered from 1 to 2340, starting at the 
top of each scene. 

(2) Samples are numbered from 1 to 3240 (approximately) starting 
at the left of each scene. 

(3) Pixels are referenced by the coordinate pair (line-number, 
sample -number) . 

1.4 EARTH- ORIENTED COORDINATES 

Earth-oriented coordinates of a control point may be in any of 
the following coordinate systems: 

(1) Geodetic Latitude and Longitude in degrees or in degrees: 
minutes ; seconds , to the nearest 0.001 degree. 

(2) Universal Transverse Mercator (UTM) Easting and Northing in 
kilometers or in meters, to the nearest 10 nveters. The 
zone number (s) must be specified. 

1.5 CONTROL NETWORK CRITERIA 

As explained previously, accurate positional measurements of a 
representative sampling of control points throughout each scene 
are needed to compute reliable transformation coefficients for 
that scene. This network of control points should meet the fol- 
lowing criteria: 

9 The control points must be distributed uniformly within the 
scene. 

® The features used for control must be stable over time. 

o The point used for control must be well-definable, both on 
the scanner data display and on recent, reliable maps. 
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• Scanner-oriented coordinates o£ control points must be precisely 
measured. 

• Earth-oriented coordinates (geographic or UTM) o£ control 
points must be precisely measured from recent, reliable maps. 

• The control network must be internally consistent, as demon- 
strated by a mathematical adjustment yielding minimal residual 
errors . 

Section 2 gives detailed procedures for laying out such a control 
network and section 3 gives detailed procedures for refining and 
validating the network. 
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2. LAYING OUT THE CONTROL NETWORK 

This section describes the initial phase o£ control network 

establishment, in which potential control points are selected 

and their approximate scanner-oriented coordinates determined. 

Specific operations in this phase are: 

(1) Outline area covered by Landsat scene on 1:1,000,000 scale 
map (preferably on Index to Topographic Map) . 

(2) Select network of potential control points for scene from 
the Landsat film image and corresponding Index to Topographic 
Maps . 

(3) Measure approximate scanner-oriented coordinates of each 
potential control point from the Landsat film image. 

Special materials required for this phase are: 

(1) Positive transparencies of Landsat 9.5-inch system-corrected 
images (scale 1:1,000,000), channel 4 (band 7). 

(2) USGS State Index to Topographic Maps (index is at scale of 

1 : 1 , 000 , 000 ). 

(3) Set of graphic devices listed in the appendix. 

2.1 LOCATING LANDSAT COVERAGE 


First, outline the area covered by the Landsat scene on a 
1:1,000,000 scale map. Within the U.S.A., the appropriate USGS 
State Index to Topographic Maps is recommended for this purpose. 
The coverage outline may be readily transferred from a 9.5-inch 
system-corrected image, either by eye or by tracing directly onto 
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the base map using a light table.* Once outlined, divide the area 
into four equal strips corresponding to the coverage of the four 
CCT's. 

Next, position and tape a i;et of the ERTS-1 MSS scales Csee fig- 
ures 3 and 4) to the system-corrected image of channel 4 (band 7) 
as follows: 

(1) Position the scale marked 18.6 CM = 3,240 SAMPLES (see fig- 
ure 3) over the image, lining up the left edge (sample 0) 
v/ith the left edge of the image, and the right edge (sample 
3,240) with the right edge of the image, and tape to the 
image . 

(2) Position the scale marked 18.1 CM = 2,256 LINES (see figure 4) 
on top, lining up line 42 with the top of the image and line 
2,298 with the bottom of the image, and tape it. 

The scales may not fit the image exactly due to the slight varia- 
tions in image dimensions. If this is the case, try to achieve 
the best average fit. When positioned properly over the image, 
the scales will appear similar to figure 5. Note that they are 
not oriented at right angles, but match the angle made by the edges 
of the image. 

The overlaid image will be used in selecting potential control ■ 
points as described in the next part of this report, as well as 
in measuring the approximate control point coordinates (para- 
graph 2.3). 


*In the event that the scene outline cannot be located precisely 
from the channel 4 (band 7) image, the channel 2 (band 5) image 
may be found more useful for this purpose. 
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2.2 SELECTING POTENTIAL CONTROL POINTS 

All potential control points must lie in areas covered by USGS 
map quadrangles of either the 7.5- or 15-minute series, as deter- 
mined from the USGS State Index to Topographic Maps.* Areas 
covered by 7.5-minute quadrangles are preferable to those covered 
by 15-minute quadrangles for the following reasons: 

(1) The 7.5-minute quadrangles are closer in scale to computer- 
generated displays of control points, and thus easier to 
visually correlate. 

(2) The 7.5-minute quadrangles are generally more recent and 
more reliable. 

(3) Geographic or UTM coordinates can be measured more easily and 
mor^T accurately from the larger-scale 7.5-minute quadrangles. 

Select the network of potential control points for each scene 
using features associated with water bodies visible on the system 
corrected image of channel 4 (band 7) (figure 1). 

Some features which may provide good control are; 

• sharp bend in a stream 

• intersection of two streams 

• major bridge across a stream 

• identifiable small lake 

• small island in a largei' lake 

• small strait 


*Many parts of the world are not co^vered by either 7.5-minute or 
15-minute quadrangle maps. However, the procedure described in 
this manual can be used with other maps; i.e., county highway 
maps. The largest scale and most accurate maps should be used. 
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Shoreline features in relatively flat terrain should not be used 
for control since they usually fluctuate with changes in water 
level. 

Six \vell-distributed control points are required for each scene 
when more than one of the four CCT's is being processed. At least 
one of these points must be near each corner of the scene and 
one near the center. In addition, each of the four CCT strips 
must contain at least one of these points. In figure 6, control 
points 1, 2, 3, 4, 5, and 6 provide an acceptable minimum, network 
for the scene. 

When processing only one CCT strip of a scene, a minimum network 
for the CCT strip may be obtained with a total of only five points 
within the strip (one at each corner and one in the center) . In 
figure 6, control points 5, 6, 7, 8, and 9 provide an acceptable 
minimum network for CCT strip number 4. 

Increasing the number of control points for each scene above 6 
will generally increase the accuracy of mapping, provided the 
points are evenly distributed. 

More potential control points should be selected than will be 
required, since several may prove to be invalid or unusable. A 
good rule of thumb is to select two potential control points for 
each valid control point desired. 

After selection of the potential control points, the next step 
is to measure their approximate coordinates using the ERTS-1 MSS 
scales overlaid on the 9. 5 -inch image. 

2.3 MEASURING APPROXIMATE CONTROL POINT COORDINATES 

Estimate the scanner-oriented coordinates of potential control 
points using the ERTS-1 MSS scales previously taped over the 
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9,5-in.ch system-corrected image as discussed in paragraph 2.1. 

For example, the picture element (pixel) in the lower right-hand 
corner of the image is located at line number 2298 and at sample 
number 3240. Following the Landsat convention of giving line 
number before sample number, we identify the scanner-oriented 
coordinates of this pixel as 2298, 3240. 

Scanner-oriented coordinates should be interpolated to the nearest 
five or 10 lines and samples. The error in these estimated coor- 
dinates W'ill generally be less than 80 lines and 50 samples. 

The following information for each potential control point should 
be recorded on a list or map: 

• control point number 

• approximate scanner- oriented coordinates 

• quadrangle name and map series 

After completion of the above, the potential control points are 
subjected to additional processing and tests before the control 
point network is considered to be established. These are 
described in the next section of this report. 
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3. EVALUATING THE CONTROL NETWORK 

This section describes the final phase of control network estab- 
lishment, in which the potential control points which were 

selected and measured in section 2 are refined and validated and 

their exact coordinates are determined. 

The specific operations in this phase are to: 

(1) produce a computer-generated display for each potential 
control point. 

(2) visually correlate the computer-generated display of each 
control point with its representation on a map. 

(3) measure exact scanner-oriented and earth-oriented coordinates 
of each control point. 

(4) adjust the control network, evaluate residual errors, and 
make necessary corrections. 

The special mateiials required for this phase are listed below: 

(1) List or map produced in the initial phase giving approximate 
scanner-oriented coordinates for all potential control points 

(2) Landsat system-corrected MSS tapes. 

(3) uses State Index to Topographic Maps (or equivalent) . 

(4) uses 7.5- or 15 -minute quadrangle Cor equivalent) for each 
potential control point. 

(5) Supply of blank Control Point Worksheets. 

(6) Set of graphic devices listed in the appendix. 


3.1 DISPLAYING POTENTIAL CONTROL POINTS 


Each Landsat pixel represents an area of approximately 1,7 acres 
on the surface of the earth. The pixels with MSS channel 4 
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(band 7) radiance values between 0 and 4 generally fall entirely 
over fresh, clear water bodies. Pixels with somewhat higher 
radiance values in channel 4 C^and T) typically fall over turbid 
water, over wetlands, or partly over water and partly over land. 

A display of pixels with channel 4 C^>and 7) radiance values from 
0 to approximately 14 will include most major lakes and streams, 
which is useful for selecting control points, but may also contain 
some features of questionable usefulness. 

Displays of all potential control points should be generated with 
the PICTAB program (see volume 2} . When executing the program in 
batch mode the window displayed for each potential control point 
should include "buffers” of at least 80 lines and SO samples 
around its estimated coordinates, to allow for errors in the esti- 
mate. A demand terminal user may achieve the same result by 
generating several displays with various spacings and origins . 
Regardless of mode, the following rules apply: 

(1) Always start with displays of channel 4 Cband 7j radiances 
between 0 and 14. 

(2) Always use displays with line and sample spacing of 1 for 
exact correlation and control point measurement. 

The computer- generated displays of the potential control points 
are next compared visually with topographic quadrangle maps of 
the areas containing the points to determine which points are 
suitable for inclusion in the control point network. 

3. 2 VISUAL CORRELATION OF CONTROL POINTS 

The user should compare the channel 4 (band 7j computer-generated 
display of each potential control point with its representation 
on a uses 7.5- or 15 -minute quadrangle map to verify that: 

(1) The feature used for control has not changed between the date 
it was mapped and the date of the Landsat overpass. 
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(2) The control point itself is well defined on both the computer- 
generated display and the map. 

This comparison is most easily performed (for maps at the scale 
of uses 7.5-minute quadrangle sheets) by overlaying the computer- 
generated display and the quadrangle map on a light table. Once 
they have been shifted and rotated to bring the control point into 
coincidence, a pointed instrument may be used to simultaneously 
prick the location on both. (This marked location should be 
exactly in the center of a print character on the display.) 

All control points should be correlated, marked, and labelled. 
Inconclusive and incorrect points should be discarded. 

The network of potential control points which have been validated 
should now be examined. These valid control points should be 
distributed evenly throughout the scene, with a minimum of one 
point near each of the four corners, and one point near the 
center (in addition, each strip must contain at least one point). 
If these conditions are not met, additional control points must 
be selected and validated. 

Figures 7 through 11 illustrate some of the situations which may 
arise in evaluating potential control points. The computer- 
generated displays in these figures were produced from CCT’s for 
the following two Landsat scenes : 

• E-1132-16535, strip 3, December 2, 1972, covering the 
Howard County, Texas, vicinity. 

• E-1073-16244 , strip 3, October 4, 1972, covering the 
Houston, Texas, vicinity. 

Each figure contains seven reduced computer-generated displays 
at the upper right, identified as T=9 through T=1S. The "T" 
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refers to the maximum channel 4 (band 7) radiance value displayed. 
A bold arrow indicates the minimum channel 4 Cband 7) radiance 
value displayed at full size on the left side of the figure. 
Channel 4 (band 7) radiance values in these full-size displays 
are symbolized as follows: 


Channel 4 Display 

Radiance Symbol 


0-9 


0-9 


10 A 

11 B 

12 C 

13 D 

14 E 

15 F 


The lower right side of each figure contains the part of a USGS 
quadrangle map corresponding to the area covered in the computer 
generated display. Note that the area covered is "skewed" due 
to the effects to earth rotation. 
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Control point evaluation (5 of 5) 
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After completing visual correlation of the potential control 
points as shown on computer- generated displays and topographic 
maps, the next step is to measure the exact earth-oriented and 
scanner-oriented coordinates of the points which satisfy the 
criteria set forth above. These exact coordinates, together with 
attitude (pitch and roll) data from the computer-generated dis- 
plays, serve as input to computer program CONTROL for mathematical 
adjustment of the control point network (part 3.4 below and 
volume 2) . 

3.3 MEASURING EXACT CONTROL POINT COORDINATES 

Exact scanner-oriented coordinates and geographic or UTM coor- 
dinates must be measured for all valid control points. A Control 
Point Worksheet (figure 12) should be used to record the measured 
coordinates for all control points. A partially completed work- 
sheet is shown in figure 13. Note that latitude and longitude are 
entered in degrees to the nearest 0.0001®. Also note that the 
worksheet contains a table to facilitate conversion of 2.5-minute 
parallels and meridians from degrees/minutes/seconds to decimal 
equivalents in degrees. These conversions are required when using 
the geographic nomoscales explained later in the paragraph. 
Scanner- oriented coordinates should be entered on the worksheet 
to the nearest line number and nearest sample number. Every valid 
control point should be plainly marked and numbered on both the 
computer-generated display and the USGS quadrangle map and 
described briefly on the worksheet. 

Each display generated by PICTAB contains line scales at the left 
and sample scales at both the top and the bottom. In addition, 
the displays usually contain internal tick marks every 10 lines 
and samples. By creasing or cutting the display just above the 
bottom sample scale, this scale can be positioned on a control 
point and used (in conjunction with the line scale) to read off 
its exact scanner coordinates. 
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Blank geographic control point worksheet. 
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Exact UTM coordinates can be measured from the quadrangle maps 
with either a metric scale or a special coordinate template, by 
first constructing the grid lines to the left and below the con- 
trol point Casing blue ticks in the map margin as a guide) and 
then measuring the incremental distances from these grid lines. 

Exact geographic coordinates can be measured from the USGS 7.5- 
and 15-minute quadrangle maps using a set of geographic nomoscales 
(see figures 15 and 16) . These nomoscales can be used to measure 
latitude and longitude anywhere in the western hemisphere between 
latitudes 23° North and 55° North.* Detailed instructions for the 
calibration and use of the nomoscales are contained on the scales. 
Figure 17 illustrates the calibration of a nomoscale and figure 18 
illustrates the use of a nomoscale to measure geographic 
coordinates . 

The exact scanner-oriented and earth-oriented coordinates for all 
potential control points accepted as meeting the criteria in 
paragraph 3.2 above are measured and entered on the Control Point 
Worksheets. The network of individually valid control points 
must now be adjusted mathematically and any residual errors deter- 
mined by use of the CONTROL program as discussed below. 

3.4 RESIDUAL ERRORS IN CONTROL POINT ADJUSTMENT 

The CONTROL program (see volume 2) corrects the network of control 
points for attitude and then uses a least-squares procedure to 
develop a set of transformation coefficients between earth- 
oriented coordinates and scanner-oriented coordinates which gives 
the best fit. This program then compares each control point's 
measured scanner-oriented coordinates with the scanner-oriented 
coordinates computed from the transformation coefficients. 

*Tfiese~ latitudes cover only the 48' conterminous states. The 
nomoscales must be redesigned for use in other latitudes or 
hemispheres. 
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Residual errors between measured scanner-oriented coordinates and 
computed scanner-oriented coordinates are printed to indicate 
how well the control points have been correlated and measured. 

Residual errors should meet the following criteria: 

• Every individual control point must have a residual error less 
than 250 meters. 

• The entire control network must have a root mean square (RMS] 
error less than 125 meters. 

In addition, the network should meet the following overall criteria 

• The network must cover more than 50% of the scene. 

• The network centroid should be within 20 kilometers of the 
scene center. 

• Control points should be evenly distributed (not "bunched up"). 

When the network contains a large number of control points (approx- 
imately 12 or more) evenly distributed in both lines and samples, 
an error in coordinates for a single control point should result 
in a larger residual error for that control point than for any 
other and therefore, that control point should be reevaluated or 
discarded. With fewer or poorly distributed control points, the 
point with the largest residual error may not necessarily be the 
control point whose coordinates are in error, and, therefore, all 
control points should be reevaluated. 

In any case, if the preceding criteria are not met, the offending 
control point (s) must be corrected or replaced and the network 
adjusted again. 

Once individual and RMS errors meet the criteria, the control 
network has been validated for use in producing precisely regis- 
tered maps from Lands at MSS digital data. 
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APPENDIX A 
GRAPHIC DEVICES 

• ERTS-1 MSS Scale (18.1 CM = 2,256 lines) 

• ERTS-1 MSS Scale (18,6 CM = 3,240 samples) 
» Geographic Nomoscale (1:24,000) 

• Geographic Nomoscale (1:62,500) 
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